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Abstract 



, The spin-dependent structure function of the deuteron, g\D, is calculated within 

a covariant framework. The off-shell structure of the bound nucleon gives corrections 
to the convolution model at a level of half a percent for x below 0.7, increasing to 
more than five percent at larger x. Overall, the dominant source of error comes from 



the lack of knowledge associated with the deuteron .D-state, which may introduce an 
uncertainty in the neutron spin structure function, g\ n , extracted from deuterium 
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There is currently much discussion about the interpretation of the results from the 
Spin Muon Collaboration Q and SLAC E142 ||| experiments on the polarized deuteron 
and helium structure functions. Of particular interest is the neutron structure function, 
gi n . A precise knowledge of this is necessary to verify the Bjorken sum rule — when 
combined with the previously measured proton structure function g lp || |], [|, In 
addition, an independent measurement of g\ n is important for determining the flavor 
singlet combination of polarized quark distributions — the first moment of which, in the 
naive parton model, is just the fraction of the spin of the nucleon carried by quarks. 

Before the nuclear data can be used for these purposes, it is essential to account for 
any nuclear effects that may arise when extracting guy from g\A- Nuclear corrections to 
polarized deuteron and helium structure functions, such as those due to Fermi motion, 
shadowing or meson-exchange currents, have been considered by several authors |7], |8|, |9|, 
10| , |TT[ [T2|j . Early attempts || to describe nuclear effects in the deuteron, based on time- 



ordered perturbation theory in the infinite momentum frame, suffered from the problem 
that the deuteron wavefunctions are not known in this frame. Subsequent analyses ||, 1(| 
all adopted the so-called convolution approach, in which the nuclear structure function is 
a one-dimensional convolution of the structure function of a free nucleon with the nucleon 
momentum distribution in the nucleus. Efforts to incorporate relativistic effects in the 



deuteron were made in Refs. [|11], |12| , however also within the confines of the convolution 
model. 

In contrast with the earlier work, in this letter we demonstrate that in a covari- 
ant treatment, inclusion of the full off-mass- shell structure of bound nucleons leads to a 
breakdown of convolution for spin-dependent nuclear structure functions. A convolution 
component can, however, be extracted from the full result by selectively taking on-shell 
limits for the bound nucleon structure function, and neglecting relativistic components of 
the nuclear wavefunction. Off-shell corrections to the spin-averaged structure function of 



the deuteron, F2D, were examined in Ref. ||13||, and found to be about 1-2% for x < 0.9 
(c.f. heavy nuclei or nuclear matter, where the off-shell effects can be considerably larger 
fl4l , |T5ll). Since the absolute values of g\D and g\ n are considerably smaller than F 2 d 
or i*2 n , we may expect nuclear corrections to be of greater relative importance for the 
spin-dependent structure functions. It is of some importance therefore that the issue of 
off-shell effects in g 1D be seriously addressed. 

In the present analysis we restrict ourselves to the valence component of polarized 
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structure functions, thereby avoiding confronting the issue of the axial anomaly fl6|j , 
which could divert attention from our main interest. Within the impulse approximation, 
deep inelastic scattering from a polarized deuteron is described as a two-step process, in 
terms of the virtual photon-nucleon interaction, parametrized by the truncated antisym- 
metric nucleon tensor G^ u (p, q), and the polarized nucleon-deuteron scattering amplitude, 
A(P,p, S). The antisymmetric part of the deuteron hadronic tensor can then be written 
as: 

M D W°(P, S, q) = i^- q e, ua/3 q a (V g 1D (P, q) + (V - fr^) 92d(P, q)\ 

d A p 



J-JLtc[A{P,S, P ) G^q)} 2tt5 ((P - pf - M 



(2tt)< 

where P, p and q are the deuteron, off-shell nucleon and photon momenta, respectively, 
and Md is the deuteron massQ The vector S a (S 2 = —1, P ■ S = 0) is defined in terms of 
deuteron polarization vectors e™ such that ]18[ S a (m) = —i e af3Xp e™* e™ P p /Md, where 
m = 0, ±1 is the spin projection. 

In analyzing nucleon off-shell effects, it will be convenient to expand the truncated 
nucleon tensor G^ u in terms of independent basis tensors: 

G^u(p,q) = q a (V (V75 G {p) + ^ 5 G {q) ) + 7^75 G {l) + ■■•), (2) 

where the coefficients Gu\ are scalar functions of p and q, and the dots (. . .) represent 
terms that do not contribute to the g\ structure function in the Bjorken limit, as well as 
those that vanish in the massless quark limit (which we use throughout). 

The polarized structure function, gm, can be extracted from Wj? v by considering the 
polarization combination (m=+l) - (m=-l). Using the fact that (e^ £3* — £3* ) = 
—i £\ P af3 P X S P /M D (since e+* = — e~), the deuteron-nucleon amplitude A can be written: 

A = ~2M^ €xpaf3 P X S P (p' — M)~ 1 T a (p) (P— p'—M) T^(p) (p' — M)' 1 , (3) 



where T a (p) is the relativistic deuteron-nucleon vertex function [T9|. In the massless 
quark limit only the pseudovector components of A are relevant: A = 7s7a-A a . In terms 
of A x and G({\, the gm structure function (per nucleon) is: 

9lD ^ = 4tt 2 MdS ■ q I dy ^ ( A ' q ( P ' q G(q) + G(7) ) 



+ p-q A-pG ip) ) , (4) 



1 For a discussion of some problems associated with using the impulse approximation for g2A see 
Ref.Q. 



where x = Q 2 /2P ■ q is the Bjorken scaling variable and y = p ■ q/P ■ q is the light-cone 
fraction of the deuteron momentum carried by the struck nucleon. 

The analogous hadronic tensor for an on-shell nucleon is obtained by tracing G^ v with 
the nucleon spin-energy projector: 



2MW^(p,q) 



Tr 



+ M) (l + 7 5 ft G^faq) 



(5) 



where s is the nucleon polarization vector, and by setting p 2 = M 2 and y — 1 in the 
coefficients Guy The structure function g^ can then be identified as gijsr(x) = giN{x,p 2 = 
M 2 ), where 



9in 



(6) 



The presence of the Gu,) term in Eq.(^j) (which does not appear in g 1N ) means that 
simple factorization of the fully relativistic gm into nuclear (*4 A ) and nucleon (Gm) parts 
is not strictly possible. This however would be required for convolution, as assumed in 
Refs. |], [K| |TTJ . Nonetheless, by writing A x in terms of relativistic deuteron wavef unctions, 

(see also [p0|), we can show that the Gu,) term is of order 

Indeed, all of the non- 



as calculated for example in 

(f/c) 3 compared with the first term in the integrand in Eq. 
factorizable corrections to convolution can be shown to be of higher orders in v/c. This 
is easily seen by separating the " + " component of A x into an "on-shell" part, which is 
proportional to the non- relativistic deuteron wavefunctions (see Eq.(|T2|) below), and an 
off-shell component, A + = Aq N + Aq FF , where (in the deuteron rest frame): 
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with p z = |p|cos0, E p = ^M' 2 + p 2 and cos6> = (yM D -po)/|p|. In Eqs.(|7|) and (|) 
u and w correspond to the S- and D-state deuteron wavef unctions, while v t and v s are 
the relativistic triplet and singlet P-state contributions. (In our numerical calculation we 
use the model of Ref.|]19[ with wavefunctions corresponding to the pseudovector ttNN 
interaction.) 

Using Eqs.(§)-@ we can then decompose gin into a convolution component plus 
off-shell corrections: 

giD(x) = C ^A/(y) g 1N (-) + 5 {N) g 1D (x) + ^g lD {x) + S^g 1D (x), (9) 
where || 

Am = J*P + *S(p)s(y-£y (10) 

can be identified with the difference of probabilities to find a nucleon in the deuteron with 
momentum fraction y and spin parallel and antiparallel to that of the deuteron. For a 
deuteron with polarization m — +1, AS(p) corresponds to the spectral function: 

AS(p) = tf +1 (p) S z * +1 (p) 5 ip -M D + E p ) , (11) 

where ^f m {p) is the usual (normalized) non-relativistic deuteron wavefunction, and S z is 
the z component of the nucleon spin operator. In terms of ^f m (p), Aq N can be written: 

A£> N = 8n 3 M D M AT tf +1 (p) S z * +1 (p), (12) 

where M = J d\p\ p 2 (u 2 + w 2 ) ensures that the normalization agrees with that of the 
relativistic calculation. The function Af(y) then satisfies Jq 1 dy Af(y) = 1 — 3/2 lud, 
where up = Jd|p| p 2 w 2 / Af is the non-relativistic D-state probability. 

We should point out that the definition of the convolution component in Eq. @ is not 
unique. For example, some authors take the non-relativistic limit in the argument of the 
energy-conserving 5-function, 5 (po — M — Ed + p 2 /2M), where Ed = — 2M is the 
deuteron binding energy. Furthermore, in order to make a comparison of our full result 
with previous calculations 0, we have included in Eq.fllOD the "flux factor" (1 +p z /M), 
which was introduced in Ref . || by analogy with unpolarized structure functions. However, 
this prescription also differs according to various authors || [10], ^2|. On the other hand, 
the full, relativistic result in Eq.(|]) is exact (within the impulse approximation) and 
contains no such ambiguities. 
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The three non-factorizable corrections in Eq.(0) are proportional to powers of |p|/M. 
The first two corrections, 
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arise from the off-shell components associated with the deuteron-nucleon vertex and the 
nucleon structure function, respectively. In the latter g°{j \x/y,p 2 ) = giN(x/y,p 2 ) — 
gm(x/y). The 5^ correction is given by: 

5 {G) gi D (x) 



M D M 



dy dp 



(M D -2£ P )| L-^ 



3M D p z 2 
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(15) 



To estimate the size of the relativistic corrections requires a model of the nucleon func- 
tions Guy While the scaling behavior of Gu) can be derived from the parton model, their 
complete evaluation requires model-dependent input for the non-perturbative parton- 
nucleon physics. For this purpose it is useful to analyze the problem in terms of relativistic 
quark-nucleon vertex functions, as described in Refs. |T3|, |i~4f . We use an ansatz in which 
a suitable set of spin S — and S = 1 vertex functions is chosen firstly to parametrize 
the unpolarized valence nucleon (viz. uy + dy and dy/uy) and deuteron {uy + dy) data. 
This fixes the normalization and momentum dependence of the vertex functions. The 
same set is then used to obtain the S = and 1 polarized valence distributions Ago and 
Aqi. A simple way to relate these to the polarized distributions Auy and Ady is via the 
SU(4) spin-flavor symmetric relations |^TJ: Auy = 3Ag /2 — Aqi/G and Ady = — Ag x /3. 
Of course our formal results do not rely on the use of SU(4) symmetry — these rela- 
tions merely provide a convenient way to parametrize the polarized quark distributions. 
Indeed, to reproduce the correct large-x behavior of the proton and neutron structure 
functions requires that the various spin-flavor quark distributions have different asymp- 
totic x-dependence [22|, p3[ , which necessarily breaks SU(4) symmetry. Within the present 
approach this is achieved by using different Dirac structures and momentum dependence 
for the S = and S = 1 vertices. 
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We find the polarized and unpolarized data can be fitted with the structures J$q(p, k) 
and /;$q(p, k) for the scalar vertex, and 7 5 7 Q $i(p, k) for the pseudovector, where k is 
the quark momentum. The momentum dependence in the vertices is parametrized by the 
simple form: $(p, k) = N(p 2 ) ■ k 2 /(k 2 — A 2 ) n , with N(p 2 ) determined through baryon 
number conservation. A best fit to the experimental nucleon distributions at Q 2 = 10 
GeV 2 (when evolved from Qq ~ (0.32 GeV) 2 using leading order evolution^]) is obtained 
for cut-off parameters = 1.0 GeV and Aq = 1.1 GeV, and exponents n£j = 2.0 and 
Hq = 2.8, for the two scalar vertices. These contribute to the total scalar distribution as: 
Aqo(x) = r Aq^{x) + (1 — r) Aq^(x), with r = 0.15. The parameters for the pseudovector 
vertex are Ai = 1.8 GeV and n\ = 3.2. The mass parameters associated with the 
intermediate spectator states are taken to be m (i) = (p — k) 2 = 0.9(1.6) GeV. 

With these parameters the first moments of the polarized valence distributions in the 
proton are Jq 1 dx Auy(x) = 0.99 and Jq 1 dx Ady(x) = —0.27, which saturates the Bjorken 
sum rule: Jq dx (Auy(x) — Ad v (x)) = qa- The x-dependence of the polarized proton 
structure function xgi p (x) = x(4Auy(x) + Ady(x))/18 is shown in Fig.l. In the valence 
quark dominated region (x > 0.3) the result agrees rather well with the SLAC, EMC and 
SMC proton data || f|, [J. A negatively polarized sea component at x < 0.3 would bring 
the curve even closer to the data points. 

For the deuteron, the structure function calculated from Eq.([|) is also shown in Fig.l 
(scaled by a factor 1/2). The agreement with the SMC data in the valence region 
is also quite good. Because the structure function is not very sensitive to explicit in- 
dependence in the quark-nucleon vertex functions, we take N(p 2 ) to be constant. The 
numerical values of these normalization constants are fixed through valence quark number 
conservation in the deuteron to be (1.2, 0.6, 0.8)% smaller (for the $q, ^o> vertices 
respectively) than in the case of the free nucleon. 

The resulting ratio, gm/giN, is displayed in Fig.2. For large x it exhibits the same 

characteristic shape as for the (unpolarized) nuclear EMC effect, namely a dip of ~ 7-8% 

at x ~ 0.6 and a steep rise due to Fermi motion for x > 0.6. For small x it stays below 

unity, where it can be reasonably well approximated by a constant depolarization factor, 

1 — 3/2 ud, as is typically done in data analyses [[[]]. Also shown in Fig.2 is the ratio of 
2 While a next-to-leading order analysis is important for a precise determination of the Q 2 -dependence 
of g\ and the Bjorken sum rule p^] , the present treatment is perfectly adequate for the purpose of 
evaluating the relative sizes of the nuclear corrections. 
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the convolution ansatz (Eqs.(|9|) - (|TT|)) to the full calculation (dashed curve). As one can 
see, this ansatz works remarkably well for most x, the only sizable deviations occurring 
for x > 0.8. 

To obtain a better idea about the precise origin of the off-shell effects, we plot in Fig. 3 
the individual off-shell corrections defined in Eqs.([l~3l)-(|l5|), relative to gin. For x < 0.8 
each of the corrections is of order 0.5% or less. As x — > 1, however, the convolution- 
breaking off-shell effects increase by more than an order of magnitude, and will need 
to be accounted for if one is to obtain accurate information on the x — > 1 behavior of 
the neutron structure function extracted from deuterium data. We have also made the 
calculation for the case of pseudoscalar ttN coupling ||19|| , where we find the total off-shell 
correction to be roughly twice as large p5| . However pseudoscalar coupling is usually 
considered to be less realistic than the pseudovector form. The corrections may also 
depend on the model for the off-shell nucleon structure function, although the fact that 
the off-shell effects become large at x ~ 1 is largely independent of the x — > 1 structure 
of the on-shell nucleon input [25 . 

Finally, as an estimate of the uncertainty, Agi n , introduced into the extracted g\ n 
through neglect of various nuclear effects, we plot in Fig. 4 the quantity Agi n (x) = 
^9id(x) (^R(x) — R(x)^j /R(x)R(x), where R(x) = gio(x)/ 9in(x) is obtained in the full 
model, while R(x) is determined through the convolution formula (solid) and via the 



ansatz R(x) = 1 — 3/2 lu d (dashed), with uj d = 4.7% [[19). Neglecting off-shell corrections 
then leads to an uncertainty |A#i n \ ~ 0.002 in the absolute value of gin-, which is still 
quite small. For comparison we also illustrate the error associated with using a differ- 
ent value for ujd in R(x) = 1 — 3/2 ud, namely 5.8% (dot-dashed) as used in the SMC 
analysis Jl|. A larger value for uou would have the effect of shifting the constant curve 
in Fig.2 down below the "full" result, thereby changing the overall sign of Ag lra in Fig.4. 
Thus, depending on the precise value of oj£> used, the procedure of applying a constant 
depolarization factor may lead to an absolute uncertainty |A#i n | ~ 0.01, or about 10% of 
the value of g\ n at x « 0.1 — 0.2. 

In summary, we have calculated the polarized deuteron structure function g^v within 
a covariant framework. Our analysis includes, for the first time, a detailed investigation 
of effects associated with the off-shell structure of bound nucleons — in addition to the 
more familiar Fermi motion and binding effects. The conventional convolution model 
can only be recovered from the full result by taking on-shell limits in the virtual nucleon 
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structure function and in the nucleon-deuteron interaction. In the end, our conclusion is 
that the major uncertainty in the extraction of g ln comes from the lack of knowledge of the 
deuteron D-state wavefunction. The off-shell corrections for x < 0.7 are small, of order 
0.5% out of a total nuclear effect of ~ 7%. At larger x the correction increases rapidly 
to > 5%, and will be relevant for higher moment analyses of g\ n . Although at present 
the nuclear effects still lie within the error bars of available data, in the upcoming, high- 
statistics SLAC E154 and HERMES experiments a careful analysis of all nuclear effects 
will be necessary. 
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FIGURE CAPTIONS. 

1. Valence component of the proton (solid) and deuteron (dashed) g± structure func- 
tions at Q 2 = 10 GeV 2 . The data are from Refs.|l], |3|, f|, ||, §]. The deuteron data are 
scaled by a factor 1/2. 

2. Ratio of deuteron and nucleon structure functions in the full model (solid), and 
with a constant depolarization factor 1 — 3/2 oj d (dotted), with ojd = 4.7% [pl|. Dashed 
curve is ratio of gio calculated via convolution and in the full model. 

3. Nucleon off-shell corrections to gio (dotted), S^gio (dashed), S^gio 
(dot-dashed) and the sum (solid), as a fraction of the total gm- 

4. Estimate of the error, Agi n (scaled by a factor 10), introduced into g\ n extracted 
from g\D by neglecting off-shell effects (solid), and by using a constant depolarization 
factor (1 - 3/2 u D ), with u D = 4.7% @ (dashed) and 5.8% (dot-dashed). To show 
the significance of the correction, Ag ln (xl0) is compared with the SLAC E142 g ln data 
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